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ABSTRACT 

We have conducted a series of numerical experiments using spherically symmetric, general relativistic, neu- 
trino radiation hydrodynamics with the code Agi'Ze-BOLTZTRAN to examine the effects of modern neutrino 
opacities on the development of supernova simulations. We test the effects of opacities by removing opacities 
or by undoing opacity improvements for individual opacities and groups of opacities. We find that improve- 
ments to electron capture (EC) on nuclei, namely EC on an ensemble of nuclei using modern nuclear structure 
models rather than the simpler independent-particle approximation (IPA) for EC on a mean nucleus, plays 
the most important role during core collapse of all tested neutrino opacities. Low-energy neutrinos emitted 
by modern nuclear EC preferentially escape during collapse without the energy downscattering on electrons 
required to enhance neutrino escape and deleptonization for the models with IPA nuclear EC. During shock 
breakout the primary influence on the emergent neutrinos arises from NIS on electrons. For the accretion 
phase, non-isoenergetic scattering on free nucleons and pair emission by e + e~ annihilation have the largest 
impact on the neutrino emission and shock evolution. Other opacities evaluated, including nucleon-nucleon 
bremsstrahlung and especially neutrino-positron scattering, have little measurable impact on neutrino emis- 
sion or shock dynamics. Modern treatments of nuclear electron capture, e + e~-annihilation pair emission, and 
non-isoenergetic scattering on electrons and free nucleons are critical elements of core-collapse simulations of 
all dimensionality. 

Subject headings: neutrinos — radiative transfer — supernovae: general 



1. INTRODUCTION 

As the core of a massive star collapses, electron capture 
(EC) on protons (free or within nuclei) reduces the electron 
fraction, Y e , and releases neutrinos that stream from the core. 
At densities near 6 x 10 11 g cm~ 3 , the mean-free path for neu- 
trinos of mean energy becomes comparable to the size of the 
core and the neutrinos become "trapped." The trapping of 
neutrinos effectively halts the deleptonization (reduction of 
Y e ) of the core as the emission of neutrinos through EC is 
balanced by reabsorption of the neutrinos on neutrons (either 
free or within nuclei). Neutrino opacity regulates the delep- 
tonization and sets the minimum core Yj, (Y e plus the neutrino 
fraction, Y v ), which controls the size of the homologous inner 
core, M s h oc Y^ 2 (Yahil 1983). When the density in the core 
exceeds nuclear density, the nuclear equation of state (EoS) 
stiffens, and a bounce shock forms at the sonic point that de- 
fines the edge of the homologous core. The expanding shock 
loses energy to neutrino emission and nuclear dissociation and 
stalls. The neutrino weak interactions (emission, absorption, 
and scattering) regulate energy loss by the prompt shock to 
neutrino radiation and thus affect the stalling of the shock. 

The revival of the stalled shock by neutrino heating is like- 
wise regulated by the neutrino opacities. As for photons in 
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stellar atmospheres, we can define "neutrinospheres," where 
the neutrinos of each species and energy effectively decou- 
ple from the matter in the stellar core. The neutrinosphere 
radii, are dependent on both flavor and energy with R Vc 
larger than TrL v because of the increased number of weak in- 
teraction channels for v e (neutral and charged current). The 
neutrinosphere radii increase with increasing neutrino energy 
due to the larger opacities for higher energy neutrinos. The 
emission temperatures at the various neutrinospheres effec- 
tively set the flux of neutrinos in each energy-species group 
and therefore the spectrum of neutrinos available for absorp- 
tion in the semi-transparent heating, or "gain," region between 
the proto-neutron star and the shock. 

Non-isoenergetic scattering (NIS) changes the energy and 
direction of the neutrinos, each of which plays an important 
role in all phases of the supernova evolution. Scattering to 
a lower energy can change the neutrino's environment from 
opaque to semi-transparent and allow a neutrino that would 
otherwise be trapped to escape, both increasing the luminosity 
of neutrinos available for absorption in the gain region above 
the core and enhancing lepton escape from the core. 

The earliest models of neutrino hea ting in core-collapse 
supernovae by Colg ate & White! (|l966) included only emis- 
sion by electron capture and redeposited half of the energy 
lost to core neutrino emission in the layers above the proto- 
neutron star, driving a powerful explosion of the outer lay- 
ers of the star. Subsequent spherically symmetric models 
included more detailed treatment of neutrino t ransport to 
revive the shock: gray diffusion (lArnettl 119661). two-flui d 
schemes ( Hillebran dt et all Il984t ICooperstein et al.l Il986), 
multi - group flux-limited diffu sion (Bruen rj 119751; lArnett 
[197^ iBowers & Wilson! [1981 iBruennl 119851; iMvraetaL 
1 19871). and Boltzmann t ransport dMezzacapp a & Bruenn 
1993a; Y amada et al1ll999l) . 

The neutrino opacities used evolved in complexity and 
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TABLE 1 

Neutrino Opacity Summary Table 



Interaction 


Base 


Alternate 


Model 


ve~ v'e~ 


Schinder & Shapiro (1982) 


None 


Base-noNES 


ve+ i/e+ 






Base-noNPS 




Reddv et al. (1998) 


Bruenn (1985) 


Base-ISnp 


up <-> i/p 








e — p «-v v e n 


Reddv et al. (1998) 


Bruenn (1985) 


Base-B85ea-np 


e+ra P e p 










Bruenn (1985) 


No Change 




yet <-> 


Bruenn (1985) 


No Change 




(A,Z) ^ y e (A,Z- 1) 


Lanaanke & Martmez-Pinedo (2000), 


Bruenn (1985). 


IPA 




Langanke et al. (2003) 


Fuller (1982) 




e + e — o i^P 


Schinder & Shapiro (1982) 


None 


Base-noEPpair 


A A o AAi/p 


Hannestad & Raffelt (1998) 


None 


Base-noBrems 



comp l eteness (flubbs & Schramml 119751; [Lamb & P ethicl 
119761: lYueh & Buchlerl Il976t iBludman & van Riper! 119781 
using the new weak intera ction theory (Glashow; 11961 



Weinberg 1967; Salam 1968) concurrently with relevant ex- 
perimental results including: the detection of ne utral-current 
intera ctions in neutrino-nucl eon scattering (lHasert et alj 
119731) and the discovery of r~ dPerl et al.ll 19751) . whose neu- 
trino partner, v T , was widely assumed to exist, but was not d e- 
te cted until the turn-of - the-m illenium ( Kodama et al. 2001). 

iSchinder & Shapirol (119821) and IBruennl (11985b assembled 
comprehensive neutrino opacity sets that included energy 
coupling in the scattering (NIS) and we r e orie nted toward 
numerical implementation. The IBruennl (Q985) opacity set 
("B85" hereafter) has been widely adapted both in content 
and form and has widely served as the "canonical" opac- 
ity set. The B85 opacity set includes emission, absorption, 
and isoenergetic scattering (IS) on heavy nuclei, a-particles, 
and free nucleons; NIS on electrons; and v D-pan emissio n 
from e + e~ annihilation. (Though derived in IBruennl (119851) . 
NIS on positrons is often omitted from simulations using B85 
opacities.) 

In the more than 25 years since the compilation of the B85 
opacity set, the search for opacity improvements (driven in 
part, previously, by the possibility that spherically symmet- 
ric models might explode if sufficiently detailed) has contin- 
ued. Some of the newer neutrino channels identified, devel- 
oped for simulations, and widely adopted include: nucleon- 
nucleon bremsstrahlung (Hannestad & Raffelt Hl998l and a re- 
lated inelastic ^-nucleon scattering channel); more kinemati- 
cally com plete i^-nucleon emission, absorption, and sca ttering 
opacities dBurrows & Sawverlll998t IReddv et al.lll998l) : the 
pair-flavor conversion process dBuras et al.l l2003): and tabu- 
lated EC ra tes using ensembles of nuclei with detailed leve l 
structures (Langanke et al] 120031: Uuodagalvis et al.l 12010b . 
Other refinements and corrections to opacities that have been 
added to many simulations inc lude: weak m agnetism for in- 
teractions with free nucle ons dHorowitzl 12002b: ion-ion cor- 
relations between nu clei (Bowers & Wilson 1982J iHorowitj 
1 1997b lltoh et al.l 120041) : and changes from the effective mass 
of the nucleons in dense matter (cf.. IReddv et al.lll999l) . For 
dense matter, from where nuclei become correlated through 
nuclear matter, the next frontier in the computation of v- 
nucleon/nucleus interactions is the development of EoS tables 
with consistent neutrino opacities for e mission, absorption , 
scattering, and neutrino-pa i r processes (IReddv et al.l 119981 : 
iMartfnez-Pinedo et al.ll2012l:lRoberts & Reddyll2012l) . 

Updated n eutrino opacity sets have been assembled by sev- 
eral authors (lBurrowsll200ll iBuras et al.l 120061 iBruenn et all 



l2006t ITentz et al.ll2012l) . Prior studies in spherical symmetry 
have examined th e impact of the addition, or modification, of 
a single opacity dMezzacappa & Bruennl[l993bt IBuras et al.1 
| 2003t iMesser et al.ll2003b iHix et al.ll2003b iMarek et al.ll2005b 
Langanke et al. 2008) and of multiple, s imultaneous , opac- 
ity changes in spherical symmetry (IBuras et al.l 120061 
i Lentz et all 120121) and in axisvmmetrv (2D) ( IBuras et al.l 
l2006t iMuller et al.ll2012l) . This is the first study detailing the 
effects of each of the opacity changes made to create a mod- 
ernized opacity set. In this paper, we start from the full set 
of opacities used in lLentz et al.l (120121) and test not only each 
component of the opacity change in that paper and each opac- 
ity upgrade relative to the B85 opacities, but also reexam- 
ine the omission of neutrino-electron scattering (NES) and 
all NIS opacities and omission of each pair-source opacity. 
In each of our simulations we retain at minimum the B85- 
formulation of scattering on nucleons and nuclei and emis- 
sion and absorption on nucleons and nuclei to ensure that the 
total opacity is not radically changed. While the differences 
found among the tests of opacity removal are generally con- 
sistent with prior single-opacity studies, we find that opacity 
changes in the context of a detailed opacity set can have differ- 
ent impacts than changing the same opacity in a less complete 
opacity set. This contextual effect is most prominent for NES 
during collapse, where the previously identified role of NES 
in enhancing neutrino escape and core deleptonization dur- 
ing collapse by downscattering neutrinos to lower energies 
is muted by detailed EC on nuclei. We identify, within our 
modern opacity set, critical opacities needed for reliable com- 
putation of the shock dynamics and neutrino emission during 
the collapse, shock breakout, and accretion phases, as well as 
opacities of little impact on the simulation or the observational 
neutrino properties. 

2. NUMERICAL METHODS AND INPUTS 

All models in this paper are computed using the par- 
allel version of Ag/Ze-BOLTZTRAN, a code for general- 
relativistic, s pherically-symmetric, ne utrino radiation hy- 
drodynamics (Liebendorfer et al. 2004) with extensions de- 
scribed here. 

2.1. Agi'/e-BOLTZTRAN 

Ag//e-BOLTZTRAN is a combination of the 
general-relativistic (GR) hydrodynamics code Agile 
dLiebendorfer et all 12002b and the neutrino tr ansport 
code BOLTZTRAN dMezzacappa & Bruennl Il993ab 
iMezzacappa & Messerl 11999b ITiebendor fer et al.l 120041) . 
Agile solves the complete GR spacetime and hydrodynamics 
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TABLE 2 
Model Summary Table 



Bounce properties Post-bounce peak 



Model 


Core mass (./Vfgh) 


central p c 


central Y e 


central Y^, 


shock radius 


u e -luminosity 




Mq 


10 14 gcm- 3 






km 


Bethe s" 1 


Base 


0.430 


3.234 


0.2448 


0.2804 


161 


408 


Base-noNIS 


0.431 


3.234 


0.2453 


0.2811 


150 


478 


Base-noNES 


0.430 


3.234 


0.2450 


0.2807 


158 


481 


Base-ISnp 


0.431 


3.233 


0.2451 


0.2808 


153 


404 


Base-noNPS 


0.430 


3.233 


0.2448 


0.2804 


160 


408 


IPA 


0.554 


3.824 


0.2843 


0.3331 


159 


432 


IPA-noNIS 


0.618 


4.239 


0.3099 


0.3712 


148 


449 


IPA-noNES 


0.608 


4.162 


0.3056 


0.3647 


159 


454 


IPA-ISnp 


0.554 


3.831 


0.2849 


0.3339 


149 


430 


IPA-noNPS 


0.551 


3.825 


0.2843 


0.3331 


159 


432 


Base-noEPpair 


0.431 


3.233 


0.2448 


0.2804 


183 


407 


Base-noBrems 


0.430 


3.216 


0.2443 


0.2798 


163 


407 


Base-noPair 


0.435 


3.216 


0.2443 


0.2798 


185 


410 


Base-B85ea-np 


0.431 


3.239 


0.2452 


0.2808 


159 


393 




0.2 0.4 0.6 0.8 1.0 
Enclosed Mass [M ] 




0.2 0.4 0.6 0.8 1.0 
Enclosed Mass [M ] 




0.2 0.4 0.6 0.8 1.0 
Enclosed Mass [M ] 



1.2 



E 




0.50 



0.2 0.4 0.6 0.8 1.0 
Enclosed Mass [M ] 



1.2 



0.45 



> 

S3 
O 



<3 0.40 



C3 



0.35 



-J 



0.30 



H 0.25 
W 




0.2 0.4 0.6 0.8 1.0 
Enclosed Mass [M ] 



1.2 




0.0 0.2 0.4 0.6 0.8 1.0 
Enclosed Mass [M ] 
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FIG . 1 . — Properties of models at core bounce, where bounce is defined as the maximum compression of the central density during the launching of the b ounce 
shock. Models shown are: Base (black; all opacities), Base-noNES (blue; without NES) with the other NIS opacity variation models discussed in Se ction|3.1.1| 
indistinguishable from models Base and Base-noNES at bounce and omitted for clarity; the NIS opacity variation models with the IPA EC from Section [3.1.2| IPA 
(orange; all NIS opacities), IPA-noNIS (red; without NIS opacities; no NES, no NPS, nucleon IS), IPA-noNES (green; wi thou t NES), but not IPA-ISnp (nucleon 
IS), whi ch is indistinguishable from model IPA at bounce and omitted for clarity. The pair opacity test models (Section |3. 3) and improved nucleon EC model 
(Section|5]2j are also indistinguishable from model Base and omitted for clarity. The panels are: radial velocity (upper left), density (upper center), entropy (upper 
right), temperature (fcT, lower left), net electron (or proton) fraction (Y e , lower center, solid lines), net lepton fraction (Yl = Y e + (n V(1 — rii/ e )/n| : , aryona> 
lower center, dashed lines), and pressure (lower right). All quantities are plotted relative to enclosed rest-mass in Mq. 
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equations implicitly in spherical symmetry on a dynamic, 
moving grid. The moving grid allows us to adequately re- 
solve whole collapsing core, including the shock, using only 
0(100) radial zones. Enhancements include the use of a TVD 
(total variation di minishing) advection sch eme in the hydro- 
dynamics solver (Liebendorfe r et al.l 120051) . which improves 
the accuracy of advection, and the use of 8m as the grid 
coordinate rather than the enclosed mass (lFischeretal.f 2010. 
Section 2.1), which improves numerical accuracy when mass 
zones are small and density gradients are large. BOLTZTRAN 
(Mezzacappa & Bruenn 1993a; Mezzacappa & Messer 1999; 
Liebendorfer et al. 2004) solves the GR neutrino Boltzmann 
equation using the method of discrete ordinates (Sn) with a 
Gauss-Legendre quadrature. Here we use an 8-point angular 
quadrature and 20 logarithmically-spaced energy groups, 
with group centers from 3 to 300 MeV. The discretization 
scheme is designed to simultaneou sly conserve lepton num - 
ber and energy as described in Liebendor fer et al.l (2004). 
Since we do not include any physics to distinguish between 
muon- and tau-fiavored leptons, we use the combined species 
*V = v r) and v^ T = {u^, v T }. 

2.2. Opacities and other inputs 

For all models we use the n uclear, electron, and pho- 
ton equations of state (EoS) of Latti mer & Swestyl (119911) 
with the bulk incompressibility of nuclear matter k s = 
220MeV[] This m atches the current ex perimental value of 
k s = 240 ± 20 MeV dShlomo et al.l2 006) better than the value 
of 180 MeV more commonly used with LS EoS in the past, 
though the value of k s has been shown to be of little con- 
sequence during the earl y phases of core-collapse supernova 
evolu t ion shown here (ISwestv et al.l Il994t iThompson et al.l 
120031: iLentzet all 120101) . Matter outside the "iron" coreQ is 
treated as an ideal gas of 28 Si that "flashes" instantaneously to 
nuclear statistical equilibrium when the temperature exceeds 
kT > 0.47 MeV. 

The stellar progenitor used for all mode ls reported here is 
the solar-metalicity progenitor of Woos ley & Hegerl 

(120071) . We have mapped the inner 1.8 Mq of the progenitor 
onto 108 mass shells of the adaptive radial grid. 

The base opacity set inc ludes emission, a bsorption, and 
scattering on free nucleons (Redd y et al.lll998]): isoenergeti c 
scattering on heavy nuclei and a-particles (Bruenn 1985); 
scatteri ng of neutrinos on electro ns (NES) and positrons 
(NPS) ( Schinder & Shapiro] 119821); production of neutrino 
pairs from e + e~ annihilation JSchinder & Shapiro! 119821) 
and n ucleon-nucleon bremsstrahlung (IHannestad & Raffelt 
1998); and electro n capture (EC) on n uclei using the 
LMSH EC table of lLanganke et al] (120031). w h ich u tilizes 
the EC rates of lLa nganke & Martfne z-Pinedol ( 2000) and 
ISampaio et al.l (120021) . The full angle and energy exchange 
for scattering between the neutrinos and electrons, positrons, 
and nucleons is included, while scattering on nuclei is isoen- 
ergetic. Bremsstrahlung and e + e~ annihilation are the only 
sources of v^ T and v^ T . 

In this paper we conduct numerical experiments where in- 
dividual opacities in the modernized base set are replaced 
with alternatives, or removed. For NES, NPS, and the e + e~- 

7 We use the latest version of the Lattimer & Swesty (1991) EoS, 
version 2.7, which is available for download from its authors at 
http : / /www .astro . suny sb .edu/dswesty/lseos. html 

8 The Fe-core is defined as the inner core that is in nuclear statistical equi- 
librium (NSE) or where the mass fraction of Fe-peak nuclei exceeds 0.5. 



annihilation and bremsstrahlung pair sources we have no al- 
ternatives so these opacities are tested by removal. For emis- 
sion, absorption, and scattering on nucleons, and EC on nu- 
clei, the base opacity com ponents a re rep laced individually 
by the simpler versions in iBruennl (119851) . The alternative 
to the LMSH EC table, which includes detailed emission 
rates over an ensemble of nuclei, is an independen t parti- 
cle approximation (IPA) dFullerlll982t IBruennl [19851) . which 
cuts off when the mean neutron number of the heavy nuclei 
N > 40. The alternative to the NIS nucleon scattering opaci- 
ties of lReddv et al.ldl 998) are the more approximate IS equiv- 
alents from IBruennl (119851) . which include the phase space, 
but not the recoil of the nucleons. The alternative to the neu- 
trino emission and abso rption on free nucleon opacities of 
iReddv et al.l (fl998) from |Bruenn| (1985) uses an approximate 
phase space factor and omits nucleon recoil effects. Ion-ion 
correlations and weak magnetism are omitted from ou r opac- 
ity set and tests. The scattering on nuclei remains the Bruenn 
( 1985) form for all models. The base opacities and their al- 
ternatives along with the models testing each alternative are 
summarized in Table Q] 

3. RESULTS 

A previous paper dLentz et al.l 12012b examined the effects 
of removing significant sections of the available modern opac- 
ity set with noticeable consequences. However, the effects of 
individual opacity changes were not isolated and the opacity 
comparison was conducted using Newtonian hydrodynamics 
and gravity and 0(v/c) transport. In this paper we start from 
the same general relativistic model with the full modern opac- 
ities and examine the effects of the various opacities in detail. 
The tests are organized into three groups: 1) non-isoenergetic 
scattering opacity tests, 2) emission/absorption opacity tests, 
and 3) pair opacity tests. Each modern opacity is removed 
or replaced with an alternative, individually, with additional 
models in which groups of opacities are removed or replaced 
with alternatives. The models are summarized in Table [2] 
The configurations at bounce of all distinguishable models 
are plotted in Figure Q] We limit our simulations to the first 
150 ms after bounce when shock radius, proto-neutron star 
radius, and thermodynamic profiles are reasonable approxi- 
mations to multi-dimensional models. 

3.1. Non-isoenergetic Scattering Comparisons 

The r e is a rich literature (lBruennlll985l iBruenn & Haxtonl 
19911; iMezzacappa & Bruennl Il993bt ISmit et al l 119961; 
Thom pson et al.l 120031) on NES effects in core-collapse 
demonstrating that, despite being a relatively small con- 
tributor to the total scattering opacity, neutrino-electron 
scattering (NES) plays an important role in determining core 
deleptonization during collapse. The dominant opacities that 
control the flow of neutrinos and lepton number from the core 
are energy dependent, with larger opacities (lower mean free 
paths) for higher neutrino energies. Energy downscattering 
by neutrinos on electrons (or other constituents) lowers 
the energy and therefore the optical depth of the scattered 
neutrino. This opens a natural channel for enhanced neutrino 
escape that is artificially suppressed when NIS opacities are 
excluded from the core-collapse model. 

We retest the effects of energy downscattering on collapse 
and the subsequent evolution of the collapsed core by re- 
moving each source of NIS individually from our reference 
model Base. We also compare to a model that includes no 
NIS opacities. Surprisingly, given the previous studies, none 
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FIG. 2. — Shoc k trajectories in km, versus time after bounce, for all models 
in Section 13.1.1 1 The models plotted are Base (black; all opacities included), 
Base-noNIS (red; without NIS opacities; no NES, no NPS, nucleon IS), Base- 
noNES (green; without NES), and Base-ISnp (blue; nucleon IS). Shock posi- 
tion is computed by bisecting the pair of mass shells with the largest negative 
radial velocity gradient —dv r /dr. 
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FIG. 3. — Comovi n°-fram e neutrino luminosities measured at 400 km for 
all models in Section l3.1.1l Colors are as in Figure [2] Electron neutrino, u e , 
luminosities are represented by solid lines, v e -luminosities by dotted lines, 
and i/^t -luminosities by dashed lines. For the heavy lepton flavor neutrinos, 
L MT = L M = L T , is the luminosity of a single species and not a sum. P^t- 
luminosities are indistinguishable from v^ T -luminosities, and omitted from 
this figure. The luminosities are in Bethe s , where 1 Bethe = 10 51 ergs. 

of these changes affected the collapse ph ase, th ough differ- 
ences do emerge after core bounce (Section fe.l.H . To sort out 
these differences, we repeated the full set of numerical exper- 
iments after removing a key opacity unavailable in the previ- 
ous studies — the LMSH electro n capture table. The results 
using the IPA EC (Section \3. 1.2b follow the general expecta- 
tions of the previous studies. The analysis of the difference 
made by the choice of EC opacity during collapse (Section 
13 - 1 . 3b provides a cautionary warning about the coupled na- 
ture of the neutrino opacity contributions to collapse and the 
viability of explosions. 

3.1.1. NIS comparisons using LMSH EC table 

For this set of tests we compare a model with our full opac- 
ity set (Base) to models without electron scattering (Base- 
noNES), without positron s cattering (Base-noN PS), replac- 
ing the nucl eon scattering o f iReddv et al.l (Q~998) with the IS 
equivalent of lBruennl(ll985l) (Base-ISnp), and to a model with 
all three of these changes (Base-noNIS). At bounce we find no 
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FIG. 4. — Comoving-frame neutrino RMS energies, (£„)rms = 
(/ dn dE E 4 F/ f dfxdE E 2 F) 1 / 2 , measured at 400 km for all models in 
Section [3. 1.1 1 RMS energy is computed over number density, not number 
flux. Colors are as in Figure [2] Line styles are as in Figure [3] 

noticeable differences among the models. (See Figure Q] and 
Table |2]) As the shock reaches its largest extent (Figure [2), 
two splittings are visible. (The model without positron scat- 
tering, Base-noNPS, is omitted from the plots and discussion 
as it is indistinguishable from the Base model.) The smaller 
splitting reaching 2-3 km by rs 80 ms post-bounce is created 
by omitting electron scattering. The larger splitting, origi- 
nating at « 30 ms post-bounce, is 10-15 km and is created 
by replacing the NIS nucleon scattering with the IS equiva- 
lent. In both cases, the larger radius is obtained by the model 
with more NIS opacity. The splitting in the shock radii is re- 
flected in the luminosities for the same epoch (Figure[3]i. The 
choice of nucleon scattering opacity makes the largest differ- 
ence, with the nucleon NIS-containing models (Base, Base- 
noNES) having larger luminosities for all neutrino species, 
«5-6 Bethe s _1 for each species at 150 ms after core bounce, 
due to smaller tota l scattering opacity fo r the nucleon N IS of 
IReddv et all (119981) than the nucleon IS o flBruennl d 19851) . per- 
mitting easier escape of trapped v e v e from the core. A much 
smaller difference in luminosities exists for the omission of 
NES. The larger v e - and ^-luminosities for models using 
the NIS nucleon scattering sustain larger shock radii through 
higher gain-region net heating rates. For the smaller splitting 
in luminosity and shock radius due to removing NES, the cor- 
relation between luminosity and shock radius is weaker, being 
inverted in the case including nucleon NIS. We were unable 
to isolate a cause for this small difference. For the two mod- 
els using nucleon IS, the model with NES (Base-ISnp) has 
slightly higher luminosity and shock radius than the model 
without NES (Base-noNIS); but, for the two models with nu- 
cleon NIS, the model with NES (Base) has lower luminosity 
and higher shock radius relative to the model without NES 
(Base-noNES). ? have proposed an analytic proportionality 

for the shock radius, R s h oc (L l/ (El)) i / 9 Rj,l '> where R vs 
is the neutrinosphere radius. If we consider only the effect of 
f e P e -luminosity we would find a relation closer to i? sn oc L v 
for the scattering models showing l arger d ifferences in lumi- 
nosities in this section and Section 13.1.21 indicating that the 
increase in neutrinosphere radii from lower opacities (for the 
more luminous models) are also playing a role in the shock ra- 
dius. There are no detectable differences among these models 
in the (E^Jrms and (-EpJrms (Figure H} during the accre- 
tion epoch. During the accretion phase, the various neutri- 
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FIG. 5. — Shoc k trajectories in km, versus time after bounce, for all models 
in Section l3~j~2l Models are IPA (orange; all NIS opacities), IPA-noNIS (red; 
without NIS opacities; no NES, no NPS, nucleon IS), IPA-noNES (green; 
without NES), and IPA-ISnp (blue; nucleon IS). 

nospheres and the entire gain region are within a gas com- 
posed primarily of electrons and free protons and neutrons, 
and the nucleon scattering opacity improvements are reflected 
in the radiation and hydrodynamic quantities. In addition 
to permitting non-isoenergetic scattering, the improved nu- 
cleon opacities account for nucleon phase blocking, recoil, 
etc., which also make modest alterations to the total scatter- 
ing opacity. 

During breakout we see a different behavior, with NES be- 
ing the more important NIS opacity. The breakout burst L„ s is 
approximately 480 Bethe s _1 for models without NES (Base- 
noNES, Base-noNIS) and approximately 405 Bethe s _1 for 
models with NES (Base, Base-ISnp). The same grouping 
applies for (-E v Jrms> which is up to 2 MeV higher during 
breakout and up to 1 MeV higher before bounce for models 
without NES. The breakout burst represents the passage of the 
shock through the neutrinospheres, and the material above the 
shock is primarily composed of heavy nuclei. Therefore, scat- 
tering on free nucleons is of less importance to the neutrino 
spectrum during breakout. The electrons above the shock 
downscatter the energies of the escaping neutrinos and lower 
the (KJrms and, therefore, the total luminosity. 

The most dramatic difference seen in these models is in the 
quantity (£ , ^ 1t )rms (Figure|U dashed lines), with significant 
differences across all the models. At 150 ms post-bounce 
the increase in (-B^ t )rms relative to the Base model is ap- 
proximately 0.5 and 1.5 MeV, respectively, for models Base- 
noNES and Base-ISnp, while the increase is nearly 7 MeV for 
model Base-noNIS, with all NIS scatterings removed. This 
demonstrates the non-linear nature of the opacity changes. 
Energy lost by via scattering thermalizes the spectra of 
v^r reducing (-E^ t )rms an d serves as an important source 
of heati ng b etween the v^ lT - and ^ e -neutrinospheres. (See 
Section 13.31 for further discussion.) In model Base-noNIS 
there are no scattering processes remaining to thermalize the 
spectrum of v^ r after emission. NIS from either nucleons or 
electrons is enough to significantly lower (-E^^rmSi with 
both required for the full effect. The lack of variation in 
{E Vc p c ) rms a ft er breakout for these models demonstrates that 
absorption followed by emission plays the role of an "effec- 
tive downscattering" in thermalizing the neutrino spectra. 

3. 1.2. NIS comparisons using IPA EC 
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FIG. 6. — Comovin g-fram e neutrino luminosities measured at 400 km for 
all models in Section [3.1.2l Colors are as in Figure |5]Line styles are as in 
FigurefS] The luminosities are in Bethe s , where 1 Bethe = 10 ergs. 
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FIG. 7. — Comovin g-fram e neutrino RMS energies measured at 400 km for 
all models in Section [3.1.2| Colors are as in Figure [5] Line styles are as in 
Figure [3] 

The lack of effects at bounce fro m reduced NIS opaci- 
ties sta nds in contrast to the results of Mezzacapp a & Bruennl 
(1993b) who used an earlier version of our code with the 
same NES opacity. Therefore, we repeated our NIS subtrac- 
tion exper iments with th e L MSH EC table replaced by the 
IPA EC o flBruennl d 19851) that lMezzacappa & Bruennl dl993bl) 
used. The replacement of the LMSH EC table with the IPA 
EC has the largest effect on the bounce configuration of any 
pair or emission/absorption opacity alternatives. (See Fig- 
ure Q] Thos e d iffer ences will be discussed individually in 
Sections [3.2l and r3~3l ) At core bounce (Figure[T|i, the primary 
differences among these models are attributed to the inclu- 
sion or omission of NES. The homologous core mass, M s h, 
shifts from 0.544 M Q for model IPA with all NIS opacities, to 
0.608 M Q for model IPA-noNES with NES removed. There 
is a slightly larger shift in M s h to 0.618 M@ for model IPA- 
noNIS without the remaining NIS opacities. (The model with- 
out NPS, IPA-noNPS, is indistinguishable from the model IPA 
throughout and is omitted from the plots and discussion.) Cor- 
related to the increase in M sn is Yj,, which increases from 
0.333 for IPA to 0.365 for IPA-noNES. There is a correspond- 
ing increase in Y e from 0.284 to 0.306 and in net neutrino 
number, Y v = Y L - Y e , from 0.048 8 to 0.0591, consistent 
with the effects of NES described by Mezzacapp a & Bruennl 
( 1993b), where scattering by NES moved neutrinos to lower 
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energies with lower opacities and an easier path to escape. 
In the models without NES, neutrinos are not scattered to 
lower energies by electrons, and are not aided in their escape. 
During collapse, the core consists primarily of heavy nuclei 
and electrons, therefore NIS on free nucleons has little effect, 
as we can see from these models. The model with nucleon 
IS and including NES, IPA-ISnp, is indistinguishable from 
model IPA at bounce. The models with higher core Y e have a 
correspondingly higher core density, pressure, and post-shock 
temperature. 

During breakout, the inclusion or omission of NES remains 
the primary difference among the models. The shocks (Fig- 
ure[5]l for these models with the IPA EC on nuclei launch more 
vigorously and from shallower depths in the gravitational well 
than the shocks for the LMSH EC models (Figure |2}. The 
vigor of these launched shocks creates a slight oscillation in 
shock radius for the models with NES (IPA, IPA-ISnp) and a 
prominent shock "ringing" in the models without NES (I PA- 
noNES, IPA-noNIS), also seen bv lThompson et al.1 (12003b for 
a model without NES. The breakout luminosity for v e peaks 
somewhat higher at L Ve ~ 450 Bethe s _1 for the models 
without NES and w 430 Bethe s _1 for the models that include 
NES. (See Table[2]) The breakout burst luminosity in v e also 
drops faster for the models without NES and exhibits oscilla- 
tions in all luminosities and (E v )-rms (Figure |7). The RMS 
energies for v e are also higher in the pre-bounce phase for 
the models without NES. During breakout the shock passes 
through the various ^ e -neutrinospheres, and in the case of the 
models without NES, the shock oscillates through the neutri- 
nospheres. 

At 20 ms after bounce the shock trajectories cross (Figure^) 
and by 40 ms after bounce the primary difference between 
models is the use of NIS or IS nucleon opacities. Like for 
the models with the LMSH table, the shock radius in the later 
ep ochs is ~ 1 0- 1 2 km larger for models using the nucleon NIS 
of lReddv et all d 19981) (IPA, IPA-noNES) than those using the 
nucleon IS of |Bruenn| (UHl (IPA-ISnp, IPA-noNIS). Within 
these pairs of models and during this epoch, the difference 
caused by NES is even smaller for the IPA EC models than 
the LMSH EC models. This difference in shock radii is also 
reflected in the luminosities (Figure [6), with larger luminosi- 
ties for models with nucleon NIS than for those with nucleon 
IS. Lower scattering opacity for nucleon NIS again enhances 
escape of trapped v e v e , increasing their luminosities, the heat- 
ing from absorption, and the shock radius. The smaller shifts 
in luminosity within these nucleon scattering pairs due to their 
differences in the inclusion or omission of NES are not re- 
flected in the shock radii. 

The late-time behavior of (£i/)rms (Figure 13 for these 
IPA EC models is similar to the LMSH EC models. After 
» 40 ms, differences in (E Vc )hms an d (£'p ( ,)rms are unde- 
tectable for these models. The peak (E„ T }rms (dashed lines) 
during breakout changes from 18 MeV for model IPA with 
all NIS opacities included, to approximately 23 MeV for the 
two models with one NIS opacity missing (IPA-noNES, IPA- 
ISnp), to 42 MeV for model IPA-noNIS with all NIS opaci- 
ties omitted. For the IPA-noNIS model, there are no energy- 
exchanging scatterings to alter the spectrum, and what we see 
is the emission spectra of all emitted neutrinos integrated over 
the semi-transparent pair-emitting region. Inclusion of either 
NES or nucleon NIS is enough to push the observed v^ T - 
spectrum most of the way toward the values in model IPA 
with all of the NIS opacities. A similar effect is also seen 



at 150 ms after bounce when the models with one missing 
NIS opacity (IPA-noNES, IPA-ISnp) have (E^Jrms that is 
< 2 MeV larger than for the IPA model, while for the model 
without any NIS opacities (IPA-noNIS), (S^Jrms is 6 MeV 
larger than for model IPA, demonstrating the non-linear inter- 
play of the NIS opacities on the thermalization of the f AtT C' AiT 
spectra. 

3.1.3. Interaction between NIS and EC during collapse 

To illustrate the interplay between nuclear electron cap- 
ture and non-isoenergetic scattering during collapse, we plot 
the occupation number (where values of 1 represent a com- 
pletely occupied phase-space) of v e as a function of energy, 
the zeroth moment of the distribution function, for three pre- 
bounce epochs. In Figure [8] we compare the spectra of mod- 
els with the LMSH EC (Base, Base-noNIS; black lines) to 
models using the IPA EC (IPA, IPA-noNIS; orange lines), 
for both cases: including the full set of NIS opacities (Base, 
IPA; solid lines) and omitting all NIS opacities (Base-noNIS, 
IPA-noNIS; dashed lines). During collapse the central den- 
sity, p c , serves as a useful "clock" for comparing different 
models, and the spectra are plotted for p c — 10 11 , 10 12 , and 
10 13 g cm -3 (top, middle, and bottom panels, respectively). 
At ,o c = 10 11 g cm -3 (upper panel) the spectra are relatively 
similar, though the IPA EC models have fewer neutrinos at 
low energies and slightly more neutrinos in the high-energy 
tail. At p c = 10 12 g cm -3 (center panel) there is a clear sep- 
aration between the models. The overall phase-space occupa- 
tion has increased for all models, though for the lowest energy 
neutrinos in model IPA-noNIS, it has only slightly increased. 
The high-energy tail for the IPA EC models (orange lines) is 
shifted w 40-50% higher in energy relative to the LMSH EC 
models (black lines), with the model lacking NIS, IPA-noNIS, 
having the largest shift. At energies less than « 15 MeV, 
model IPA-noNIS has fewer neutrinos than the other models, 
approximately 30 times fewer for the lowest-energy neutri- 
nos. Unlike the other models, which have fairly flat neutrino 
spectra up to the high-energy roll-off, model IPA-noNIS has 
a peak at 20 MeV. At p c = 10 13 g cm -3 (lower panel) the 
trend continues with model IPA-noNIS as the outlier. The 
high-energy tail of the v e -spectrum for IPA EC models (or- 
ange lines) is again shifted w 40-50% in energy relative to 
the LMSH EC models (black lines) with a smaller shift for 
models without NIS (dashed lines) relative to those with NIS 
(solid lines). The neutrino phase-space is (nearly) completely 
filled for E < 20 MeV for the LMSH EC models (black lines) 
and for E < 30 MeV for model IPA, while model IPA-noNIS 
reaches a peak at 30 MeV and decreases for lower energies. 
The spectra for the IPA EC models (Figure [8] orange lines) 
are consistent with previously rep orted collapse-phase spec- 
tra (Em^IilillsmtSililla). 

In the independent particle approximation, nuclear EC is 
completely shut-off when the mean neutron number N > 40, 
which occurs at densities above « 2 x 10 10 g cm~ 3 , consti- 
tuting much of the core collapse. Therefore, in the IPA EC 
models (IPA, IPA-noNIS) all of the v e emission in regions 
with density exceeding w 2 x 10 10 g cm~ 3 arises from elec- 
tron capture on free protons, which are relatively rare, and 
results in a slower overall rate of v e emission and core delep- 
tonization. (Pair emission is corres pond ingly low during col- 
lapse as we shall discuss in Section [331 ) The LMSH EC im- 
plementation emits neutrinos with a lower mean energy than 
the capture on protons, which dominates the total capture rate 
in the IPA models at densities above k 2 x 10 gem -3 , 
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FIG. 8. — Spectral neutrino occupation number (zeroth angular moment of 
the invariant neutrino distribution function) versus comoving-frame neutrino 
energ y at mass shell M = 0.25 Mq for the models discussed in Section 
13.1.31 The models are Base (black, solid), Base-noNIS (black, dashed), IPA 
(orange, solid), and IPA-noNIS (orange, dashed). The three panels are for 
the models at central densities, p c , of 1 X 10 11 g cm -3 (upper panel), 1 X 
10 12 g cm -3 (middle panel), and 1 X 10 13 g cm -3 (lower panel). 

and fills the lo w energy spectra whic h is underpopulated by 
EC on protons (Langanke et ai1 l2003l) . which we can see by 
comparing the spectral evolution of the models without NIS 
opacities for the LMSH EC (Base-noNIS) and IPA EC (IPA- 
noNIS). With IPA EC, the low-energy spectrum can only be 
filled by energy downscattering by NIS opacities, which effi- 
ciently fill the low-energy spectrum for model IPA. Using the 
LMSH EC fills the low-energy v e spectrum directly without 
energy downscattering (model Base-noNIS), so we do not see 
an enhancement in neutrino escape through the less opaque 
neutrino "window" at lower energies when NIS opacities are 
included, like we do for the IPA EC models. 

Scattering on nuclei is the dominant contributor to total 
opacity during collapse and identical in all of our models. 
Once the low-energy phase space is filled, deleptonization 
is controlled by the total opacity, and thus there are no dis- 
cernible differences in the dele ptoniz ation among the LMSH 
EC models presented in Sec tion l3 . 1.11 For the IPA EC models 
presented in Section 13.1.21 changes to the filling of the low- 
energy phase space, do affect the net escape after trapping 



begins. The differences between the models with the full set 
of NIS opacities (Base, IPA) a re a consequenc e of differences 
in EC before trapping occurs dHix et al.ll2003l 

After considering the effects of compression due to varia- 
tion in central density on the temperature and the change in 
entropy profiles due changes in M s h among the models plot- 
ted in Figure Q] there is likely an additional small increase 
in temperature and entropy inside the shock for models that 
include IPA EC and NES (orange lines). This small thermal 
increase is likely due to energy gained by the fluid from the 
scattering of v e by electrons as can be seen in difference be- 
tween the model IPA and model IPA-noNIS spectra in Fig- 
ure[8] The lack of similar effects on the entropy and tempera- 
ture by changes in scattering on models using the LMSH EC, 
suggests (as does the lack of impact on the spectra during col- 
lapse) that the NIS has very little impact on v e during collapse 
when utilizing the more complete LMSH EC implementation. 

3.2. Emission and absorption comparisons 

To test the two sources of v e (y e ) emission via electron 
(positron) capture and the inverse neutrino absorption in the 
moder n opac it y set of model Base, we replaced them with 
their iBruennl (|T985) equivalents for EC on nuclei (model 
IPA) and EC on free nucleons (model Base-B85ea-np). At 
bounce (Figure[T|i, the change in nucleon EC produces no dis- 
cernible difference between models Base and Base-B85ea-np, 
but changing the nuclear EC produces the largest differences 
at bounce of any single opa c ity rep lacement tested in this pa- 
per. As noted by iHix et al.l (120031) . the IPA EC (model IPA) 
results in less deleptonization, as the IPA turns off completely 
where the heavy nuclei have neutron numbers N > 40, of the 
inner core relative to using the LMSH EC (model Base) lead- 
ing to higher Y e , Yl, density, pressure, and temperature inside 
the bounce shock. The IPA EC also results in higher delep- 
tonization and stronger collapse outside the bounce shock, 
which can be seen in the lower Y e , higher density, and higher 
infall velocities for model IPA outside the bounce shock (Fig- 
ure [TJ where the stronger IPA for low density does not shut 
off. This will increase the ram pressure as the shock passes 
through the Fe-core. The shock forms at 0.554 Af© in model 
IPA and 0.430 Mq in model Base, and the shallower shock 
launch results in a more vigorous launch of the shock for 
model IPA (Figure |9). Approximately 70 ms after bounce 
the shock tr ajectories of t he two nuclear EC models cross (as 
reported bv lHix et al.ll2003l) due to higher ram pressure from 
higher density and infall velocities in the outer core of model 
IPA induced by stronger deleptonization. The shock trajecto- 
ries of models Base and IPA cross again late in our simula- 
tions, with the net effect that the shock position of model IPA 
is flatter during the epoch when we should expect multidimen- 
sional effects to become important. In the model using the 
B85 nucleon EC (Base-B85ea-np), the shock radius trails that 
in the Base model slightly, with the deficit reaching 5 km by 
the end of our simulations at 150 ms post-bounce. This deficit 
in shock radius for model Base-B85ea-np is reflected in an 
w 2Bethe s _1 lower u e - and ^-luminosity during the accre- 
tion phase relative to model Base, with a small (15 Bethe s _1 ) 
decrease in the breakout ^-luminosity. The shallower bounce 
shock of model IPA also results in a breakout ^-luminosity 
burst that peaks sooner and higher (by 25 Bethe s _1 ) than 
the LMSH EC model (Base), and then drops sooner creating 
a narrower breakout peak. Model IPA shows a higher v^ T - 
luminosity from bounce to approximately 80 ms post-bounce. 
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FIG. 9. — Sh ock trajectories in km, versus time after bounce, for all models 
in Section l3~2l Models are Base (black; all opacities), IPA (orange; IPA EC), 
and Base-B85ea-np (green; B85 EC on free nucleons). 
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FIG. 1 1. — Comoving- fram e neutrino RMS energies measured at 400 km 
for all models in Section 1331 Colors are as in Figure|9] Line styles are as in 
Figure [3] 
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FIG. 10. — Comovi ng-fr ame neutrino luminosities measured at 400 km for 
all models in Section 13.21 Colors are as in Figure [9] Line styles are as in 
FigurefS] The luminosities are in Bethe s" 1 , where 1 Bethe = 10 51 ergs. 

With no discernible differences in (^„ t )rms outside the nar- 
row peak during breakout (Figure ITU , the higher luminosity 
of model Base-B85ea-np implies it is emitting more v^v^- 
pairs than the other models. The (£v)rms f° r all models is 
essentially identical for each neutrino species during the ac- 
cretion epoch, with sharper and narrower peaks for model IPA 
during breakout, owing to the shallower and rapidly moving 
shock. Before bounce, model IPA shows a 1-2 MeV higher 
(£i/ 6 )rms. reflecting the shift of the hig h-energy tail seen in 
Figure|8]and discussed in Section [3.1.3l 

3.3. Pair opacity comparisons 

To test the effect of pair opacities on the supernova models, 
we compare the Base model to models omitting e + e~ annihi- 
lation (Base-noEPpair), bremsstrahlung (Base-noBrems), and 
both (Base-noPair) pair sources. Without any pair sources 
model Base-noPair also lacks v^v^. At bounce there are 
no distinguishable differences among these models, which is 
consistent with the lack of thermal e + e~ pairs to annihilate 
and free nucleons for bremsstrahlung. The shock trajecto- 
ries for the models with missing pair opacities begin to de- 
viate from model Base approximately 40 ms after bounce 
(Figure [121 . with the model Base-noEPpair shock extend- 
ing to 183 km, the model Base-noBrems shock extending 
to 163 km, and the model Base-noPair shock extending to 
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FIG. 12. — Shock trajectories in km, versus time after bounce, for the 
pair opacity test models in Section 13.31 The models are Base (black; all 
opacities), Base-noEPpair (red; without e+e — -annihilation), Base-noBrems 
(green; without Bremsstrahlung), and Base-noPair (blue; without pair opaci- 
ties). 



185 km. These are the only models in this paper that exceed 
the maximum shock extent of the Base model at 161 km. (See 
Table |2]) Comparison of these models shows that e + e _ an- 
nihilation is much more important to the shock propagation 
than bremsstrahlung. 

Unlike the previous comparisons, where accretion-phase 
shock differences are correlated with v e - and ^-luminosities, 
the luminosities (Figure [T3l for these models are anti- 
correlated with the shock radii. The differences in the 
(-£-„) rms for all neutrino species (Figure \14\ are also anti- 
correlated, with the shock with smaller differences relative 
to model Base for model Base-noBrems and larger differ- 
ences for model Base-noEPpair. In model Base-noPair the 
only source of v e is electron capture on protons and of v e is 
positron capture o n neu trons. As in the case of the EC-NIS 
interplay (Sec tion fe. 1.3b . the non-linearity of pair opacity ef- 
fects on the emission of v e v e is clear. For both luminosity and 
RMS energy the removal of bremsstrahlung had only a minor 
impact when e + e _ -annihilation was present. However, when 
e + e~-annihilation was absent, the removal of bremsstrahlung 
had a significant impact o n both. 

Previous studies dThompson et ail 120001 [2003 ; 
iBurrows et all l2000t lKeiletal.1 l2003t [Burasetall 12003 , 
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FIG. 13. — Comovi ng-fr ame neutrino luminositie s m easured at 400 km for 
all models in Section 13.31 Colors are as in Figure [12] Line styles are as in 
Figure[3] The luminosities are in Bethe s _1 , where 1 Bethe = 10 51 ergs. The 
dashed line for model Base-noPair L„ flT is omitted as this model does not 
include u UT P UT . 
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FIG. 14. — Comoving -fram e neutrino RMS energ ies measured at 400 km 
for all models in Section l331 Colors are as in Figure[l2] Line styles are as in 
Figure [5] 



2006) have stressed t he importance of b remsstrahlung as a 
neutrino-pair source. iBuras et ail (120061) plotted (see their 
Figures 21-22) each opacity for all neutrino species at two 
energies for two post-bounce epochs within our study range 
and found that bremsstrahlung was the dominant pair-source 
in the proto-neutron star, but e + e~ annihilation became 
dominant starting somewhere outside the neutrinosphere in 
each case. The v^v^ emission rates for our models are con- 
sistent with that finding, being higher in the inner core for the 
bremsstrahlung-including models (Base, Base-noEPpair) than 
the model without bre msstrahlung (Base-noB r ems). Station- 
ary transport studies (iThompson et al.l 120001 iBurrows et al. 
I2000t iKeil et alJl2003h have found bremsstrahlung to be more 
important than e + e~ annihilation to the emitted flux and 
thermalization rate of v^v^. Stationary transport studies 
drive toward a global equilibrium that is never reached during 
the dynamic shock revival phase and can overemphasize the 
region below the neutrinospheres where bremsstrahlung is 
dominant. A more relevant comp arison is with the dynamic 
models of Thompso n et al.l (120031) computed with and with- 
out bremsstrahlung like our models , though with different 
progenitor, base opacities, code, etc. IThompson et al.l (120031) 
find that the total luminosity for all heavy lepton neutrinos, 
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FIG. 15. — Net heating rates in ergs g _1 s _1 for models Base (black) and 
Base-noEPpair (red) at 30 ms after bounce. Line styles indicate net heat- 
ing by: emission/absorption by v e (solid) and P e (dotted), by f Mr P MT pair 
sources (dashed) and by NIS of v^tD^t (dash-dotted). Net pair heating 
for v & v e is approximately one order of magnitude or more smaller than the 
f/ir^/jT counterpart and is omitted from this plot for clarity. 

L x = 2L^ T + 2Lfif, increases by 8-10 Bethe s _1 starting 
from about 30 ms after bounce onward. Our models are con- 
sistent with this change, showing an increase of 2 Bethe s _1 
for each of the four heavy lepton (anti)neutrino species during 
the same epoch. The proportionate increas e due to a dding 
brems strahlung does appear larger for the IThompso n et al.l 
(2003) models as their model without bremsstrahlung has 
roughly half the neutrino luminosity for each species as 
our Base-noBrems model. Much of the difference can be 
attributed to our use of nucleon NIS and GR with the rest due 
to other differences between the opacities, progenitors, and 
codes. At later epochs, beyond the effective use of spherical 
models, the diffusion of core v^v^ to the neutrinospheres 
may enhance the relative contribution of bremsstrahlung to 
the luminosity. 

Removing pair sources reduces cooling by v^v^ emis- 
sion during the accretion phase. Figure [l5]shows the net heat- 
ing rate for models Base and Base-noEPpair at 30 ms after 
bounce, with the net cooling for v^v^ emission shown as 
dashed lines. For the Base model, the cooling by v^v^ 
emission dominates the radiative heat budget between 25 and 
40 km, with about half of the cooling returned back to the 
local thermal pool by NIS on v^v^. (In both of these mod- 
els the full set of NIS opacities are available for thermaliza- 
tion during scattering.) When e + e~-annihilation is removed 
(Base-noEPpair, red lines), both the pair cooling and related 
NIS heating are reduced and shifted deeper into the core. 
The lower effective cooling by ^ Air i' Mr -pairs in model Base- 
noEPpair results in more thermal energy being available to 
support the shock and a larger shock radius. Emission of v e v e - 
pairs is suppressed by the largely filled v e phase space. 

4. CONCLUSIONS 

In this paper we have systematically examined the effects 
of each of the updated opacities in our modernized opacity 
set over the initial 150 ms post-bounce, spherically symmet- 
ric phase of core-collapse supernovae. To summarize our pri- 
mary findings: (1) During collapse, electron capture on heavy 
nuclei dominates the emission of neutrinos and the delep- 
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tonization of the core. If modern EC rate s are used with a 
detaile d NSE composition, as in LMSH or Juodagalv is et alj 
(2010), the direct emission of low-energy neutrinos obviates 
the need for NES to fill the low-energy portion of the neutrino 
spectrum. (2) Omitting NES results in an sa 15% increase in 
the breakout-burst v e -luminosity. (3) Changes in shock for- 
mation due to deleptonization are the primary opacity-driven 
source of variation in early shock evolution. (4) During the 
accretion phase, nucleon NIS enhances the z/ e ^luminosities, 
net heating, and drives the shock further out, enhancing the 
potential of shock revival by multidimensional effects. (5) 
Cooling by 2/ A(T z/ AlT -pair emission from e + e~-annihilation re- 
moves energy from the system that could otherwise be used 
to revive the shock. (6) All of the NIS opacities (except scat- 
tering on positrons) and e + e _ -annihilation affect the neutrino 
luminosities and/or (-E f )rms during accretion. (7) Positron 
scattering shows no impact on the outcome or observables 
during our simulations. 

We have identified non-linear behaviors in the interplay 
among opacities, which illustrate that the context provided 
by the included opacities is important in evaluating individ- 
ual opacities. Some examples of neutrino opacity interplay 
include: 

• Emission from nuclear EC and energy downscattering 
by NES compete to fill the lowest energy bins during 
collapse. The escape of low-energy v e increases core 
deleptonization. The low-energy spectrum of neutri- 
nos emitted by the LMSH EC table fills the low-energy 
phase space adequately without NES; thus, we do not 
see an impact on deleptonization when NES is omitted 
as we do in the case of models using the IPA EC, which 
does not fill this part of the spectrum directly. 

• Thermalization of v^v^ by individual NIS opacities 
is not simply additive. Removing either NES or nu- 
cleon NIS results in a modest increase in (£v )rmS> 
while removing both results in a much larger increase, 
as thermalization by scattering above the emission re- 
gion is absent. 

• Neutrino emission by pair sources also exhibits satu- 
ration effects. In models including e + e~ annihilation, 
bremsstrahlung has only a minor impact on the emit- 
ted neutrino properties, but when bremsstrahlung is the 
only pair source, its removal has a much larger impact. 

We can identify from our tested set necessary neutrino- 
matter interactions required for modern supernova modeling 
in any dimension. 



• Modern nuclear EC (LMSH, I Juodagalvis et all 120101 
or equivalent) should be considered an essential ingre- 
dient in any reali stic supernova si mulation as was pre- 
viously noted by lHix et al.l ((2003). Relying on the IPA 
EC artificially alters the electron capture, deleptoniza- 
tion, and the impact of other opacities. 

• Nucleon NIS extends the shock radius via an increase 
in v e v e -luminosity. The related enhancements to cap- 
ture on free nucleons, though relatively modest in ef- 
fect, should be included for physical consistency of the 
nucleon opacities. 

• NES significantly reduces v e emission during breakout 
and contributes to thermalizing the v^v^r spectra. 

• ^r^T-pair emission by e + e~-annihilation is an im- 
portant source of cooling during the accretion phase, 
while bremsst rahlung plays only a small role (as 
also seen in Thomps on et al.l [2003) unless e + e~- 
annihilation is omitted. Bremsstrahlung may become 
more important at later epochs as trapped v^ T v^ T dif- 
fuse out. 

Omitting any of these opacities would alter the observable 
neutrino properties, and thus would introduce unnecessary 
systematic errors in the analysis of observed supernova neu- 
trino signals. 

The modern opacities discussed in this pape r are physically 
well-motivated improvements to the reference Bruenn (1985) 
opacity set. Including these improved opacities increases the 
physical fidelity of the neutrino-matter interactions in super- 
nova simulations, while omitting them risks potential system- 
atic errors in the dynamical and observational properties of 
simulated supernovae. 
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